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The breakdown of statistical homogeneity and isotropy of cosmic perturbations is a generic feature 
of ultra large scale structure of the cosmos, in particular, of non trivial cosmic topology. The statisti- 
cal isotropy (SI) of the Cosmic Microwave Background temperature fluctuations (CMB anisotropy) 
is sensitive to this breakdown on the largest scales comparable to, and even beyond the cosmic hori- 
zon. We study a set of measures, Kg (I = 1, 2, 3, . . .) which for non-zero values indicate and quantify 
statistical isotropy violations in a CMB map. The main goal here is to interpret the Kg spectrum 
and relate it to characteristic patterns in the correlation function of CMB anisotropy arising from 
cosmic topology. We numerically compute the predicted Kg spectrum for CMB anisotropy in flat 
torus universe models. The essential features are captured in the leading order approximation to the 
correlation function where Kg can be calculated analytically. The Kg spectrum is shown to reflect the 
number, importance and relative orientation of principal directions in the CMB correlation dictated 
by the shape of the Dirichlet domain (DD) of the compact space and its size relative to cosmic 
horizon. Hence, besides detecting cosmic topology, Kg can discriminate between different topology 
of the universe complementing ongoing search for cosmic topology in CMB anisotropy data. 

PACS numbers: 98.70Vc,04.20,Gz,98.80.Cq 
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In standard cosmology, the Cosmic Microwave Back- 
ground (CMB) anisotropy is expected to be statistically 
isotropic, i.e., statistical expectation values of the tem- 
perature fluctuations AT(q) are preserved under rota- 
tions of the sky. In particular, the angular correlation 
function C(q, q') = (AT(g)AT(g')) is rotationally in- 
variant for Gaussian fields. In spherical harmonic space, 
where AT(q) = J2i m a imYi m (q) this translates to a diag- 



onal 



,) = Ci8w8 mm i where Cj is the widely used 



angular power spectrum of CMB anisotropy. 

It is important to determine whether the CMB sky is a 
realization of a statistically isotropic process, or not from 
the observations themselves. We study a set of measures 
Kg (l = 1,2,3, . . .) that measure violation of statistical 
isotropy 

The detection of statistical isotropy (SI) violations can 
have exciting and far-reaching implication for cosmology. 
The realization that the universe with the same local ge- 
ometry has many different choices of global topology has 
been a theoretical curiosity as old as modern cosmology. 
Motivations for cosmic topology and their consequences 
have been extensively studied |2| . CMB anisotropy mea- 
surements have brought cosmic topology from the realm 
of theoretical possibility to within the grasp of observa- 
tions 0- 0| ■ A generic consequence of cosmic topology is 
the breaking of statistical isotropy in characteristic pat- 
terns determined by the photon geodesic structure of the 
manifold. Global isotropy of space is violated in all multi 
connected models (except S 3 jZ^). In cosmology, the 
Dirichlet domain (DD) constructed around the observer 
represents the universe as 'seen' by the observer. The 
SI breakdown is apparent in the principal axes present 
in the shape of the DD constructed with the observer 
located at the basepoint. 

In this paper we compute and study the Kg spectrum 



of SI violation arising in flat (Euclidean) simple torus 
models with a cubic, cuboidal and more generally, par- 
allelepiped (squeezed) fundamental domain. The CMB 
anisotropy in torus spaces has been well studied 0, |]| . 
We can relate the Kg spectrum to the principal directions 
normal to pair of faces of the DD, their relative orienta- 
tion and the relative importance given by the distance to 
the faces along them. Along the most dominant axes, the 
distance is minimum, and equals the inradius, R < , the 
radii of largest sphere fully enclosed within the DD [j| . 

The set of measures Kg of statistical isotropy violation 
is defined as 
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dK X g(n)C(Tlq,Kq') 
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where C(lZq, TZq') is the two point correlation between 
TZq and TZq' obtained by rotatingjy and q' by an el- 
ement TZ of the rotation group Ml. The measures 
Kg involve angular average of the correlation weighed 
by the characteristic function of the rotation group 
Xt(K) = Ejlj^mmC^) where D mm> are the Wigner 



D-functions [6|. When SI holds C(Hq, TZq') = C{q, q') 
is invariant under rotation, and eq. Q gives Kg = kq 5go 
due to the orthonormality of Xi(TZ). Hence, non-zero Kg 
for £ > measure violation of statistical isotropy. 

The measure Kg has a clear interpretation in harmonic 
space. The two point correlation C(q, q') can be ex- 
panded in terms of the orthonormal set of bipolar spher- 
ical harmonics as 

C(q, q>) = AffiYm ® Y v (q')}g M , (2) 

U'lM 

where A\jf are the coefficients of the expansion. These 
coefficients are related to 'angular momentum' sum over 
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FIG. 1: The pattern in CMB correlation in a multi-connected universe, M = M u /F is related to the distribution of 'images' of 
the sphere of last scattering (SLS) on the universal cover. The three figures illustrate this for three values of a — 45°, 60°, 75° 
in a squeezed torus. The solid parallelepiped is the fundamental domain and the dashed show its images tessellate,^ 11 . The 
solid circle is the SLS and dashed/dotted ones are its images. In each case, the radius of SLS is equal to the inradius, 7?< of 
Ai. The shaded polygon is the Dirichlet domain (DD). The closest SLS images which determine the DD are dashed (others are 
dotted). Note the hexagonal shape of DD when a 7^ n/2 and is equal sided for a = tt/3. For smaller a, the DD approximates 
an elongated cuboid. 



the covariances (ai m a*, m ,) as 

A \]f = ( a l™ a *'m>) tfml'-m' . ( 3 ) 

mm' 

where £f^/ m / are Clebsch-Gordan coefficients |jj . When 
SI holds (ai m a*, m ,) = Ci6w8 mm >, implying A e u ¥ = 
(-l) l Ci(2l + l) 1 / 2 S w S eo Smo- Aff represent the statisti- 
cally isotropic part of a general correlation function. The 
bipolar functions transform just like ordinary spherical 
harmonic function Ylm under rotation Substitut- 
ing the expansion eq. (J2J into eq. (JJJ we can show that 
Kg = J2wm \ Aw J \ 2 is positive semidefinite and is also 
given by 

Kt = 2 8 ^ 2 1 fd1lXe(Tl) ( a lrna*Vm')( a lmai, m ,) n , (4) 

lml'm f 

where (. . .) n is computed in a frame rotated by TZ 0. 

The compact spaces with Euclidean geometry (zero 
curvature) have been completely classified. In three di- 
mensions, there are known to be six possible topologies 
that lead to orientable spaces 0, 0- The simple flat 
torus, M = T 3 , is obtained by identifying the univer- 
sal cover M. u = £ 3 under a discrete group of translations 
along three non-degenerate axes, si, S2, S3: s, — > Si+nLi, 
where Li is the identification length of the torus along 
Si and n is a vector with integer components. In the 
most general form, the fundamental domain (FD) is a 
parallelepiped defined by three sides Li and the three 
angles a, between the axes (We call it squeezed torus). 
If Si are orthogonal then one gets cuboidal FD, which 
for equal Li reduces to the cubic torus. The cuboid and 
squeezed spaces which can be obtained by a linear co- 
ordinate transformation C on cubic torus [f| can have 
distinctly different global symmetry |Tl|. 

We restrict attention to the case where CMB 
anisotropy arises entirely at the sphere of last scattering 



(SLS) of radius i?» (nearly equal to observable horizon). 
Invoking method of images, the CMB correlation pattern 
on the SLS is known to be dictated by the distribution 
of nearest 'images' of the SLS on the universal cover 
The correlations are distorted even when the SLS and 
its images do not intersect (i?* < i?<). When SLS in- 
tersects its images the CMB sky is multiply imaged in 
characteristic correlation pattern of pairs of circles 
Fig. illustrates the role of the nearest SLS images and 
related DD in defining the principal directions in the cor- 
relation for a few cases. 

We compute C(q, q') for CMB anisotropy in torus 
space using regularized method of images 3]. We can 
compute the ki in real space using eq. Q or from A\f, 1 
using eq. in harmonic space. Fig. plots the Kg spec- 
trum for a number of cubic, cuboidal and squeezed torus 
spaces. We find the following interesting results : 

i. Kg — for odd £ for all torus models. This does 
not hold for compact space of non-zero curvature, e.g., 
compact hyperbolic spaces. 

ii. For cubic torus K2 = 0. K2 is non-zero for cuboidal 
and squeezed torus. This is a clear signature of non- 
cubic torus where the DD differs from the FD and has 
more than three principle axes. 

iii. For equal-sided squeezed torus, K4 , decreases as 
a decreases from 90° to 60° as i?< increases. For a < 
60° sharply increases with decreasing a as i?< decreases 
sharply yjj. 

iv. «2 = increases monotonically as a decreases from 
90°. The trend is well fit by eq. (Q , 

v. The peak of Kg shifts to larger £ for small spaces. 
The results can be understood using the leading order 
terms of the correlation function in a torus where Kg can 
be calculated analytically. For brevity, we outline the 
steps for the cubic torus. The results for cuboid and 
squeezed case is readily obtained using the transforma- 
tion C. Here we list only the results leaving details of the 



3 



2r 

1.5 ■ 
1 ■ 

0.5 ■ 
■ 



4 - 
2 ■ 


1 

0.8 r 
0.6 
0.4 ■ 
0.2 
■ 



Cubic 
(big) 



■■ R=0.4L 

I I R=0.5L 
I I R=0.6L 



Cuboid 
eql. vol. 



Ja 



Squeezed 
eql. vol. 



□ W 



■■ 45 Deg 

I I 60 Deg 
I I 75 Deg 



50 
40 
30 
20 



10 

- 
9 1 

-1 4 



9 1 



1.5 

1 

0.5 





Cubic 
(small) 



Inn i In 



R=1.0L 
R=1.5L 
R=2.0L 



10 11 12 13 



Cuboid 
R=R 



Squeezed 
R=R 



L 2= 5L 1 

=2 L 



■■ 45 Deg 

I I 60 Deg 
I I 75 Deg 



FIG. 2: The Ke spectra for flat tori models are plotted. The top row panels are for cubic tori spaces. The left panel shows 
spaces of volume, Vm, larger than the volume V* contained in the sphere of last scattering (SLS) with Vm/V* = 3.7, 1.9, 1.1, 
respectively. The right panel shows small spaces with Vm/V» = 0.24,0.07,0.03, respectively. Note that K2 = for cubic tori. 
The middle panels consider cuboid tori with 1 : 5 and 1 : 2 ratio of identification lengths. The bottom panels show ki for 
equal-sided squeezed tori with a = 45° , 60° and 75° . In the middle and bottom rows, the right panels show the case when radius 
of SLS, R, — R< the inradius of the space. Here, the SLS just touches its nearest images (see Fig.0 which is at the threshold 
where CMB anisotropy is multiply imaged for larger i?*. The cases in the left panels of lower two rows have Vm/V* = 1 and 
are at the divide between large and small spaces. 



calculation to a more comprehensive publication [8j- 

The spatial correlation function of gravitational poten- 
tial in the periodic box implied by the T 3 topology with 
cubic fundamental domain is 



£*(x,x') = L- 3 ^P$(fc n ) e 



-t2js-(x-s 



(5) 



where n is 3-tuple of integers, fc 2 = (27r/L) 2 (n ■ n), and 
the term with n • n = is excluded from the summation. 
For Naive Sachs- Wolfe CMB anisotropy arising at the 
SLS ( AT = 1/3$), the correlation function, C(q,q r ) is 
given by the spatial correlation of $ at points on the SLS 
along the two directions q and q' as 

C(q,q') = L~ 3 Y,P*(kn) e^*"'^' ( 6 ) 
ii 

where the small parameter e q < 1 is the physical distance 
to the SLS along q in units of L/2 (more generally, L/2 



where L = [LxL^L-; 



Contribution to C(g, <f) from 



large wavenumbers, |n|e » 1 is expected to be approx- 
imately statistically isotropic. The SI violation is found 
in the low wavenumbers. 



When the SLS is contained with the fundamental do- 
main around the observer, i.e., the CMB anisotropy is 
not multiply imaged on the sky, e = 2i?»/ L is a constant. 
When e is a small constant, the leading order terms in 
the correlation function eq. I© can be readily obtained 
in power series expansion in powers of e. For the lowest 
wavenumbers n| 2 — 1 in a cuboid FD torus 



C(q, q')&2j2 P^IU) cosine^Aq, 



(J) 



where Aqt are the components of Ag = q — q' along the 
three axes of the torus and Pi = L/Li. Only even powers 
of eAqi are present in the expansion eq. JJJ. This holds 
for the terms from higher wavenumbers (|n| 2 = 2 or 3) 
and explains the strictly zero ki for odd i [8|- 

In a cubic (equal sided) torus, up to the leading order 
SI violating term, the correlation is 



C(q,q')^C 



l-e 2 \Aq\ 2 + 3e^( A ^) 4 



We retain the e 4 term since the term at e 2 is explicitly 



(8) 
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rotationally invariant, hence does not contribute to the 
violation of SI. The non-zero Kt can be analytically com- 
puted to be 
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(9) 



The first non zero Kg in cubic (equal-sided) torus occurs 
at I = 4. k 4 ~ e 8 falls off rapidly as e — » for large 
spaces. 

For the cuboidal (unequal-sided) torus, the correlation 
violates SI at order e 2 



C(q,q')^C 
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The non zero Kg corresponding to the correlation eq. 1101) 
are 



observed non-zero Kg and discriminate between different 
topologies of the universe. When CMB anisotropy is mul- 
tiply imaged, the Kg. spectrum corresponds to a correla- 
tion pattern of matched pairs of circles Q. Kg have an 
advantage of being insensitive to the overall orientation 
of the correlation features. Also Kg are sensitive to SI 
violation even when CMB is not multiply imaged. The 
Afi, signature, which was not discussed here, contains 
more details of the SI violation || . 

We study a new measure of SI violation. A compan- 
ion paper describes the estimation of Ke from a CMB 
map [ij. Before ascribing the detected breakdown of 
statistical anisotropy to cosmological or astrophysical 
effects, one must carefully account for and model out 
other mundane sources of statistical anisotropy in real 
data, such as, incomplete and non-uniform sky coverage, 
beam anisotropy, foreground residuals and statistically 
anisotropic noise. These observational artifacts will be 
discussed in future publications. 
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where (3 2 = £i=iA ? - The K2 ~ e 4 signal falls off slower 
than the K4 signal in cubic space for large spaces (e — > 0). 

For the squeezed torus, the linear transformation can 
be applied to obtain the expansion up to the leading order 
SI violating term. For simplicity, we restrict to equal 
sided squeezed torus with one non-orthogonal pair of axes 
(si • S2 = cos a). The leading order approximation to 
correlation C(q, q') has the form 
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The corresponding Kg spectrum is 
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The expression for K2 explains point (iv) regarding Kg. 
Interestingly, a residual K2 remains even for very large 
space (e — * 0). 

Preferred directions and statistically anisotropic CMB 
anisotropy have been discussed in literature 10] . We 
compute SI violation of the CMB anisotropy due to cos- 
mic topology quantified in terms of the recently proposed 
Kg spectrum Q . We interpret generic features of K£ spec- 
trum arising from the shape of the Dirichlet domain as 
well as through analytic results for leading order approx- 
imation to the correlation. These results help interpret 
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For cubic torus the Dirichlet domain (DD) matches the 
fundamental domain (FD). However, for torus spaces 
with cuboid and parallelepiped FD, the corresponding 
DD is very different, e.g., hexagonal prism (see Fig0. 



